Transforming carbon nanotubes by silylation: An ab initio study 
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We use ab initio density functional calculations to study the chemical functionalization of single- 
wall carbon nanotubes and graphene monolayers by silyl (S1H3) radicals and hydrogen. We find that 
silyl radicals form strong covalent bonds with graphene and nanotube walls, causing local structural 
relaxations that enhance the sp 3 character of these graphitic nanostructures. Silylation transforms 
all carbon nanotubes into semiconductors, independent of their chirality. Calculated vibrational 
spectra suggest that specific frequency shifts can be used as a signature of successful silylation. 
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Carbon nanotubes [i. have evolved into one of the 
most intensively studied one-dimensional nanostructures. 
They owe their success in the nanotcchnology field to a 
unique combination of atomic-scale perfection, structural 
stability, beneficial thermal and electrical properties [lj. 
Their usefulness as molecular nanowires is hindered by 
two main limitations. First, reproducible formation of 
contacts with a well-defined geometry and electronic 
properties is difficult, since most materials do not wet 
the surface of the inert nanotube, but rather bond to 
ill-defined unterminated ends, causing large variations in 
the contact resistance [2J. Even more important is the fact 
that the conductance of carbon nanotubes depends sen- 
sitively on the precise way a graphene layer is rolled up 
into the tubular structure, identified by the chiral index 
(n, m) j3jlH[5]. Since it is currently impossible to produce 
nanotubes with a unique chiral index, even samples with 
a narrow diameter range necessarily contain a mixture of 
metallic and semiconducting nanotubes, reducing their 
usefulness in quantum devices. 

Here we propose that the recently [6 achieved chemi- 
cal functionalization of carbon nanotubes by silyl groups 
bears promise for the nanotube walls to form strong, 
well-defined bonds to a surrounding matrix, and for the 
majority of nanotubes to convert into semiconducting 
nanowires. Using ab initio calculations, we demonstrate 
that silyl (-S1H3) radicals form strong covalent bonds 
with the surface of nanotubes and graphene. Since sily- 
lation introduces a partial sp 3 character into the sp 2 
graphitic network, it opens up the fundamental band 
gap and converts all carbon nanotubes to semiconduc- 
tors. We show that successful silylation of nanotubes can 
be verified experimentally by monitoring specific changes 
in the vibrational spectra [6]. 

The selection of silyl radicals as functional groups has 
been motivated by the success of silane chemistry in es- 
tablishing strong, covalent bonds between polymers and 
inorganic matter [7]. The rich silane chemistry [5] is likely 
to provide ways to modify silyl radicals bonded to the 
nanotube wall in order to optimize covalent bonding of 
nanotube walls to many materials, including silicon. To 
our knowledge, there has been no prior theoretical study 



of the silyl-nanotube interaction and the effect of silyla- 
tion on the nanotube properties. 

In order to determine, whether silyl radicals may co- 
vale ntly attach to nanotubes and thereby modify their 
conductance properties in a desirable way, we used ab 
initio calculations to determine the equilibrium geome- 
try, total energy and electronic structure of silylated nan- 
otubes. Our calculations are based on the density func- 
tional theory (DFT) within the local density approxima- 
tion (LDA). We used the Perdew-Zunger [9] parameter- 
ized exchange-correlation functional, as implemented in 
the SIESTA code [ID] , and a double-^ polarized basis lo- 
calized at the atomic sites. The valence electrons were de- 
scribed by norm-conserving Troullier-Martins pseudopo- 
tentials in the Kleinman-Bylander factorized form 

In our calculations, we considered the products of the 
dissociative adsorption of silane, SifD— >SiH3(ad)+H(ad), 
yielding silyl radicals and hydrogen atoms chemisorbed 
on graphene and carbon nanotubes at different cover- 
ages. We obtained results for carbon nanotubes of differ- 
ent diameters, with chiral indices ranging from (4, 4) to 
(18,0). To describe isolated nanotubes while using peri- 
odic boundary conditions, we arranged them on a tetrag- 
onal lattice with a large inter-wire separation of 22 A. 
Depending on the SiH 3 coverage and the chiral index, we 
used supercells containing up to 4 primitive unit cells in 
the tube direction [13] and sampled the Brillouin zone of 
these ID structures by at least 7 k points. We limited 
the range of the localized orbitals in such a way that the 
energy shift caused by their spatial confinement was no 
more than 50 meV[I2]. The charge density and the po- 
tentials were determined on a real-space grid with a mesh 
cutoff energy of 200 Ry, which was sufficient to achieve a 
total energy convergence of better than 0.1 meV/atom. 

When calculating the vibrational spectrum of pristine 
and silylated nanotubes, we first optimized the reference 
structure to a high precision, with residual forces not ex- 
ceeding 0.04 eV/A. The high mesh cutoff energy value 
guarantees that the forces are translationally invariant. 
To obtain a good estimate of the Hessian matrix in the 
harmonic limit, we displaced atoms in the positive and 
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FIG. 1: (Color online) Bonding of SiH3 and H to graphitic car- 
bon, (a) Equilibrium adsorption geometry on graphene in side 
view (top panel) and top view (bottom panel), (b) Charge 
density difference Ap(r) = p(SiH3 + H/graphene) — p(SiHa) — 
p(H) — p(graphene), displayed in the plane indicated by the 
dashed line in (a), (c) Adsorption energy i? at j(SiH3+H) on 
carbon nanotubes as a function of the inverse tube diameter 
1/d. 



negative direction along each axis by the small amount 
of 0.01 A and determined the forces analytically. With 
these numerical precautions and limitations of the LDA 
force field, we found our frequencies to be dependable 
with the accuracy of down to few inverse centimeters. 
We also found our LDA- value [15] for the G-band fre- 
quency of 1630 cm -1 to agree much better with the ex- 
perimental value [Tfi] of 1580 cm -1 than previous LDA 
calculations [TT]. 

We first considered co-adsorption of SiH3 and H on 
graphene, a good model for a wide-diameter carbon nan- 
otube. Our calculations suggest that both silyl and 
hydrogen preferentially adsorb on-top of carbon atoms, 
where they interact with the initially unfilled C2p z or- 
bital. We found a weak mutual attraction between silyl 
and hydrogen, mediated by the graphene substrate. The 
different adsorption geometries in neighboring adsorption 
sites are indicated schematically in Fig. [TJa). For an 
isolated silyl-hydrogen pair, we distinguished the ortho 
geometry with the adsorbates adjacent from the para ar- 
rangement with the adsorbates on opposite corners of a 
single Cg ring. In the meta arrangement, the adsorbates 
arc neither adjacent nor opposite to each other on a single 
C6 ring. Defining the adsorption energy i? a d(SiH 3 -|-H) as 
the energy to adsorb SiH3 and H from the vacuum, we 
find E a d — —4.25 eV for ortho, —3.57 eV for meta, and 
—4.31 eV for para geometry. In view of the fact that 
the adsorption energy of a hydrogen atom on graphene is 
close to —1.4 eV|18j. we conclude that the bond strength 



between the silyl radical and graphene is unusually large, 
exceeding 2 eV. 

We find that the SiH4(ad)^SirI 3 (ad)+H(ad) dissocia- 
tion on graphene, associated with a local pyramidaliza- 
tion of the initially planar graphene near the adsorption 
site, occurs with an energy barrier of sa 1.5 eV in the para 
arrangement. The pyramidalization seen in the optimum 
adsorption geometry, shown in side view in Fig. [TJa), in- 
dicates a local change of the bonding character from sp 2 
to sp 3 , so that initially inert planar graphene is trans- 
formed into a reactive substrate capable of forming cova- 
lent bonds with adsorbates. The nature of the adsorption 
bond can be best seen by inspecting the charge density 
difference Ap, representing the charge redistribution in 
the system upon adsorption. The spatial distribution of 
Ap(r) for silyl and hydrogen adsorbed on graphene in 
the para arrangement is shown in Fig. [ljb) in a plane 
normal to the substrate. We observe no charge flow be- 
tween adsorbate and substrate, but rather a moderate 
charge accumulation in the Si-C and H-C bond regions, 
indicating the formation of covalent bonds. 

In contrast to graphene, the finite surface curvature 
of pristine nanotubes introduces partial sp 3 character, 
making these systems more reactive. This is best illus- 
trated by plotting the adsorption energy i? a( j(SiH3+H) 
as a function of the nanotube diameter in Fig. [ljc) for 
one SiH 3 and one H per supercell[T3]. As the sp 3 charac- 
ter increases with decreasing nanotube diameter, there is 
a decreasing energy cost to distort the substrate in order 
to optimize the adsorption bond, thus making adsorption 
more exothermic. Also at high coverages, we expect the 
adsorption bonds to be stronger on the average than at 
low coverages. In contrast to the first adsorbates in the 
supercell that introduce local substrate distortion, bond- 
ing of additional adsorbates benefits energetically from 
the previously enhanced sp 3 bonding character. Depend- 
ing on the adsorption arrangement within the supercell, 
we have found a strengthening of the average adsorption 
bond by up to 0.4 eV when increasing the coverage from 
one to three SiH3 and H units per supercell. 

Our results for the most stable adsorption arrangement 
on an (n, m) nanotube are given by solid data points, and 
less stable arrangements by the crosses in Fig. [ljc). With 
respect to our graphene results, the adsorption energy 
differences between the ortho, meta and para arrange- 
ments in (n, to) nanotubes may change by up to a few 
tenths of an eV per supercell, thus changing the most fa- 
vored adsorption arrangement. At nonzero coverages, we 
find that silyls and hydrogens prefer to form lines in axial 
direction, thus minimizing the strain energy in the sub- 
strate. Among the systems investigated here, we found 
silyl and hydrogen to bond more strongly to zigzag than 
to armchair nanotubes due to differences in the pyrami- 
dalization strain. 

Having established that silyl radicals and hydrogen 
form strong covalent bonds with nanotubes, we next in- 
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FIG. 2: (Color online) Electronic structure of (a) pristine and 
(b) silylated carbon nanotubes. The left panels display the 
effect of silylation on the band structure of the (6, 6) nan- 
otube. The right panels show the fundamental band gap E g 
in nanotubes as a function of the inverse tube diameter 1/d. 
Data for metallic nanotubes are given by solid symbols, those 
for semiconducting nanotubes by open symbols. 



vestigate the effect of chcmisorption on the electronic 
structure of carbon nanotubes. As mentioned earlier, 
the fundamental band gap of pristine carbon nanotubes 
varies strongly with the chiral index (n,m)[3] HJ In 
density functional calculations, the band gap value is 
known to be underestimated significantly. Consequently, 
the range of band gap values in nanotubes should be even 
larger than what is suggested by the DFT results for pris- 
tine nanotubes in the right panel of Fig. [2ja) . In a sample 
with a random distribution of chiral indices (n, m) , we ex- 
pect one third of nanotubes to be metallic jH HJ |5]. For 
better visual impact, metallic nanotubes are represented 
by solid data points and semiconducting nanotubes by 
open data points in the right panel of Fig. [2ja). The 
band structure E(k) of the pristine (6, 6) armchair car- 
bon nanotube is shown in the left panel of Fig. [2ja) as a 
representative for other metallic nanotubes. 

The effect of silylation on the electronic structure of 
carbon nanotubes is shown in Fig. (2^b) . Comparing the 
calculated band gaps in the silylated and pristine nan- 
otubes in the right panels of Fig. [2] we conclude that 
silylation has converted all nanotubes to semiconductors 
and narrowed down the range of their fundamental band 
gaps. Details of silylation-related changes in the elec- 
tronic structure can be best seen when comparing the 
electronic structure of the pristine and the silylated (6, 6) 
nanotube in the left panels of Fig. [2] Whereas the metal- 
lic character of the pristine system stems from the cross- 
ing of two linear bands at the Fermi level, enhancing 
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FIG. 3: (Color online) Effect of silylation on the vibrations of 
the (6, 6) carbon nanotube. (a) Low-frequency range of the 
vibrational spectrum, (b) Schematic of the radial breathing 
mode (RBM) in the pristine and silylated nanotube. 



the sp 3 character by silylation opens up a direct band 
gap in this system. Only in systems, where the number 
of H and S1H3 adsorbates per unit cell is not the same, 
a non-bonding low-dispersion band appears in the band 
gap near Ep. This impurity band is spin polarized and 
localized near the adsorbates. 

Obviously, the band gap of a given nanotube will 
change with changing sp 3 character [19], which depends 
on the adsorption geometry and the adsorbate coverage. 
Our band gap results in the right panel of Fig. [2^b) are 
for a low coverage of one SiH 3 and H per supercell[T3l 
of the pristine nanotube. Since the sp 3 character should 
increase with rising silyl coverage, and since DFT un- 
derestimates the fundamental band gap, our results in 
the right panel of Fig. |2jb) are to be considered a lower 
bound on band gaps in silylated nanotubes. 

Successful silylation of nanotubes can be verified ex- 
perimentally by inspecting the vibrational spectra as a 
specific signature of silyl-nanotube bonds. Reported Ra- 
man spectra of silylated carbon nanotubes [6] suggest that 
the dominant effect of silylation is the modification of 
mode intensities. No frequency shifts specifically re- 
lated to silylation have been observed, suggesting that 
the effect of silylation may be comparable to that of 
bundling In the following, we determine the vibra- 
tional spectra of pristine and silylated nanotubes in order 
to identify a spectroscopic signature of the silyl-nanotube 
bond. 

We have studied the vibrational spectra of (6,6) and 
(10, 0) carbon nanotubes containing up to three silyls and 
hydrogens per supercell[13 j . Upon silylation, we first 
expect the emergence of new vibrational modes in the 
nanotube, associated with Si-C and H-C stretch motion. 
Silylation furthermore modifies the vibrational modes of 
the pristine nanotube due to adsorbate-related deforma- 
tions and local stress, as well as the attachment of heavy 
SiH3 ligands. 

Comparing the vibrational spectra of (6, 6) and (10, 0) 
nanotubes, we find the general trends for silylation- 
related changes to be chirality independent, and fre- 
quency shifts to be generally less than 20 cm -1 . The 



4 



vibrational spectra of a pristine and a silylated (6, 6) car- 
bon nanotube, with three silyls and hydrogens in the or- 
tho arrangement forming a ring within the 9.9 A long 
supercell|13j. are compared to each other in Fig. |3ja). 
In the low-frequency part of the spectrum, silylated nan- 
otubes exhibit a mode at ~80 cm -1 , caused by the wag- 
ging motion of the chemisorbed radicals. Other modes 
induced by silylation include the high-frequency H-Si 
stretch in SiH 3 at 2000 cm" 1 and the H-C stretch of 
the adsorbed hydrogens at 2800 cm -1 . 

A vibration mode unique to carbon nanotubes is the 
Ai g radial breathing mode (RBM), shown schematically 
in the top panel of Fig. [3jb). As seen in Fig. [3]ja), this 
mode occurs at 276 cm" in the pristine (6, 6) nanotube. 
Even though the RBM mode is not well defined in partly 
silylated nanotubes, a closely related mode, shown in the 
bottom panel of Fig.[3jb), occurs at the higher frequency 
of 289 cm -1 in the silylated nanotube. Similar RBM up- 
shifts have been reported in pristine tubes and traced 
back to hydrostatic pressure [20j [21] or stress induced by 
inter-tube interaction [35] • We find an up-shift of RBM- 
like modes at lower coverages, when the axial distance be- 
tween adsorbed silyl radicals is large, and a down-shift at 
high coverages, when most carbon atoms are covalently 
bonded to a silyl. Much of this trend can be understood 
and reproduced by considering carbon atoms connected 
to silyls as virtual 'heavy' 43 C 'isotopes' within the pris- 
tine nanotube. In the low-coverage regime considered in 
Fig. |3jb), the heavy carbons move radially out of phase 
with the other carbons, causing a frequency up-shift. At 
high coverages, all carbons are subject to the usual radial 
breathing motion with the frequency down-shifted by the 
isotope effect. 

The high-frequency tangential mode, also known as 
the G band, is down-shifted in the silylated nanotube 
to 1619 cm -1 from its initial value of 1638 cm -1 in the 
pristine (6,6) nanotube. This frequency shift is almost 
independent of the silyl coverage. 

In conclusion, we used ab initio density functional 
calculations to study the chemical functionalization of 
single-wall carbon nanotubes and graphene monolayers 
by silyl radicals and hydrogen. We found that silyl radi- 
cals form strong covalent bonds with graphene and nan- 
otube walls, causing local structural relaxations that en- 
hance the sp 3 character of these graphitic nanostructures. 
With the help of silane chemistry, nanotube walls may 
form strong, well-defined bonds to a surrounding matrix, 
facilitating the transition to carbon-silicon hybrid elec- 
tronics. We found that silylation transforms all carbon 
nanotubes into semiconductors, independent of their chi- 
rality. Our calculated vibrational spectra suggest that 
specific frequency shifts can be used as a signature of 
successful silylation. 
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